Behavior of physical observables in the vicinity 
of the QCD critical end point 

Pedro Costa 

Centra de Fisica Teorica, Departamento de Ffsica, Universidade, P3004-516 Coimbra, Portugal 

o : 

Abstract. Using the SU(3) Nambu-Jona-Lasinio (NJL) model, we study the chiral phase transition 
' at finite T and Hb- Special attention is given to the QCD critical end point (CEP): the study of 

, physical quantities, as the pressure, the entropy, the baryon number susceptibility and the specific 

■ heat near the CEP, will provide complementary information concerning the order of the phase 
pLn ! transition. We also analyze the information provided by the study of the critical exponents around 
^ ! the CEP. 
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The existence of the CEP in QCD was suggested in the end of the eighties, and its 
properties have been studied since then (for a general review see Refs. [1,2]). The most 
recent lattice results for the study of dynamical QCD with Nf = 2+1 staggered quarks 
of physical masses indicate the location of the CEP at T CEP = 162 ± 2MeV, \x CEP = 
360±40MeV [3], however its exact location is not yet known once the location of the 
CEP depends strongly of the mass of the strange quark. At the CEP the phase transition 
Ph| is of second order, belonging to the three-dimensional Ising universality class, and this 

kind of phase transitions are characterized by long-wavelength fluctuations of the order 
parameter. 

As pointed out in [4, 5], the critical region around the CEP is not pointlike but has a 
very rich structure. The vicinity of the CEP is a privileged region to study the influence 
of different type phase transitions in the physical observables, namely, the pressure, the 
entropy, the baryon number susceptibility, Xb, and the specific heat, Cy. 

We perform our calculations in the framework of the three-flavor NJL model, includ- 
ing the determinantal 't Hooft interaction that breaks the Ua ( 1) symmetry, which has the 
following Lagrangian: 



Jgf = q{id. Y -m)q+ 8 ^Y,\W a V? + mY5WqY 



a=0 



+ 8D 



det [4(1+75)4] +det[q{l-y 5 )q] 



(1) 



By using a standard hadronization procedure, an effective meson action is obtained, lead- 
ing to the gap equations for the constituent quark masses from which several observables 
are calculated (we follow the methodology presented in detail in Refs. [6, 7]). 




FIGURE 1. The pressure (left panel) and the entropy (right panel) as functions of the temperature and 
the baryonic chemical potential. 



The fundamental relation is provided by the baryonic thermodynamic potential 

n(m,T)=E-TS- £ (2) 

i=u,d,s 

from which the relevant observables as the pressure P, the entropy S and the particle 
number /V, can be calculated as usually (the expressions are given in Ref. [6]). 

The nature of the chiral phase transition in NJL type models at finite T and/ or \Xg 
has been discussed by different authors [6]. For zero temperature the transition is of first 
order. As the temperature increases, the first order transition line persists up to the CEP. 
In the CEP the chiral transition becomes of second order. At higher temperatures there 
is a smooth crossover. 

In the vicinity of the CEP, the nature of the chiral phase transition influence strongly 
the behaviour of the physical observables. In the first order phase transition region 
(T < T CEF ', \Xb > Hg EP ) both, pressure and entropy, show a discontinuity which end 
at the CEP. In the crossover region (T > T CEF , jig < jlg EP ), the discontinuities of P 
and S vanish, and both quantities change gradually in a continuous way (see Fig. 1). 
The same situation can be seen in the behaviour of Xb'- in the first order phase transition 
region Xb has a discontinuity (left panel of Fig. 2). At the CEP, the phase transition is 
of second order, and the slope of the baryonic density tends to infinity which implies 
a diverging %g. In the crossover region, Xb changes gradually in a continuous way. A 
similar behavior is found for the specific heat for three different chemical potentials 
around the CEP, as we can observe from the right panel of Fig. 2. 

Focusing our attention on the critical behavior of the baryon number susceptibility 
in the vicinity of the CEP, we verify that Xb diverges with a certain critical exponent. 
Considering a path parallel to the fig-axis in the (T, /ig)-plane from lower fig towards the 
critical flg EP = 318.4 MeV at fixed temperature T CEP = 67.7 MeV, one can calculate 
the critical exponent £ using the linear logarithmic fit \\\Xb = — eln |jLig — fig EP \ + const, 
where the term const is independent of fig. The result that we obtain is e = 0.67 ±0.01, 
which is consistent with the mean field theory prediction: £ = 2/3. Since there is no 




FIGURE 2. Left panel: baryon number susceptibility as a function of fis for different temperatures 
around the CEP: T CEP = 67.7 MeV and T = T CEP ± 10 MeV. Right panel: specific heat as a function of 
T for different values of jj,g around the CEP: y% EP = 318.4 MeV and ji B = £i£ £/> ± 10 MeV. 



reason why the critical exponent should be equal for both regions, below and above 
flg EP , we also study Xb from higher fi B towards the critical Hg EP . Using again a 
logarithmic fit, the result is e' = 0.68 ±0.01 which is very near the value of £. This 
means that the size of the region we observe is approximately the same independently of 
the direction we choose in the path parallel to the /15-axis. Complementary information 
is also obtained from the study of the critical exponent of Cy [8]. 

Summarizing our discussion, we have analyzed the vicinity of the QCD critical end 
point in the SU(3) NJL model. We conclude that, in this region, the physical observables 
are strongly influenced by the nature of the phase transition. Around the CEP we 
have studied the baryon number susceptibility which is related with event-by-event 
fluctuations of /i in heavy-ion collisions. The study of the specific heat is also important 
once it is related with event-by-event fluctuations of T in heavy-ion collisions [9]. We 
also conclude that the critical exponents of Xb obtained in our model are consistent with 
the mean field value £ ~ e' ~ 2/3. 
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